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We present the full next-to-next-to-leading order (NNLO) corrections to the coefficient function for the polarized 
cross section d Aa/d Q of the Drell-Yan process. We study the effect of these corrections on the process p + p — > 
i + P + 'X' at an CM. energy yfS = 200 GeV. All QCD partonic subprocesses have been included provided 
the lepton pair is created by a virtual photon, which is a valid approximation for a lepton pair invariant mass 
Q < 50 GeV. For this reaction the dominant subprocess is given by q + q — > 7* + 'X' and its higher order 
corrections so that it provides us with an excellent tool to measure the polarized sea-quark densities. 



1. Introduction 

Deep inelastic electroproduction is very useful 
to extract information about the polarized va- 
lence parton densities Au v and Ad v . However 
almost no information exists about the gluon den- 
sity Ag and the sea-quark densities Au s , Ad s 
and As (including the anti-sea-quark densities). 
One of the processes proposed to measure the lat- 
ter densities is the Drell-Yan process or massive 
lepton pair production p + p — ► + l X' . In 
this reaction the dominant subprocess is given by 
valence-quark sea-quark annihilation into the lep- 
ton pair which continues to hold if we include 
higher order QCD corrections. Lepton pair pro- 
duction is given by the process 

p(P 1 ,S 1 )+p(P 2 ,S 2 )^ 1 *(q) +'X'. 

I 

->Z+(Ji)+Z-(Z 2 ) 
S=(Pi+P 2 ) 2 , Q 2 = q 2 = (h+h) 2 , (1) 

In the frame work of the parton model this pro- 
cess is described by the Drell-Yan mechanism 
where the protons are longitudinally polarized. 
At >/S = 200 GeV the photon dominates which 




Figure 1. Massive lepton pair production p+p 
7* + l X'. 



*Deutschos Eloktronen-Synchroton DESY, Platanenallee 
6, 15738 Zeuthen. 



2 




plotted in the range 2 < Q < 30 GcV at [i = Q 
in LO, NLO and NNLO. 



is characteristic of the RHIC energies at BNL. Let 
us first look at the transverse momentum p T and 
rapidity y distributions. They are given by 



d 3 Ac 



4na 2 d 2 AW 



dQ 2 dpx dy 3 Nc Q 2 dpr dy ' 



(2) 



where Nc denotes the number of colours. In the 
QCD improved parton model we have 



(PAW 
dp T dy 



E 



dx 2 A/ Ql (xi,^ 2 



a 1 ,a 2 =q,g 



2 ^ d 2 AW aia2 

xAf a2 (x2,^)— , 3 

dp T dy 

where the factorization/renormalization scale is 
li 2 and Af ai is the polarized parton density. 
The NLO corrections to the partonic distribu- 
tion d 2 AW ai a 2 1 dpr dy have been completely cal- 
culated in [1] (for the non-singlet part see also 
[2]). For pt > Q/2 the 0(a s ) quark-gluon sub- 
process dominates over all other reactions pro- 
vided there is a substantial gluon density. There- 
fore a measurement of the differential pr distri- 
bution is sensitive to the polarized gluon density. 
This is revealed by the the longitudinal asymme- 
try All = Aa /a where the set with the largest 
gluon density leads to the largest asymmetry. 



2. Computational details 

A different picture is shown by the invariant 
mass squared distribution given by 



dAcr 



dQ 2 3N C Q 2 S 



The structure function AW equals 



E 

a 1 ,a 2 =q,g 



/ dxi j 

J r Jt/xi 



dx 2 



xA/ ai (xiV)A/a 2 (x 2 , M 2 ) 



X A aia2 [ T, 



(5) 



Here the zeroth order quark-anti-quark process 
dominates so that this process is an excellent 
tool to measure the sea-quark density in proton- 
proton reactions. The NLO corrections to the co- 
efficient function A aia2 were calculated in [3] and 
the NNLO corrections have been recently com- 
puted in [4] . In higher order the quark-anti-quark 
process also dominates. If M is the amplitude of 
a parton-parton subprocess then the qq is given 

by 

q(Pi,Si) + q(p 2 , s 2 ) -> 7* +' X' , 
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Figure 4. The polarized quantity N(/j,//j,o) plot- 
ted in the range 0.4 < n/no < 2 with Q = 5 GeV 
and (ito = Q- 
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= -Tr ( 75 / 2 (jfe- m)M 75 A 

<?i(9i)Oi,si) + .9(^2) -> 7* 
IAA/J 2 ' 



(6) 
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X ± m) M"t 

<?l(?l)(Pl,Sl) +92(g2)(P2, 
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(7) 



s 2 ) -> 7* +' X 1 



X (^1 zb to) Af ^ 

9(Pi) +.9(P2)^7* 



(8) 



|AAf ffS | — £^2^2^20-2 



PS* 



P2 ' *2 

!_Tr ( M MlM * M" 1 "'^ 

Pi-h \ ' 



,(9) 



with the constraints 



$i ' Pi — j &i ' &i 



-1 



k-k = Q. 



(10) 



All these reactions contain the 75-matrix and/or 
the Levi-Civita tensor. If we use n-dimensional 
regularization we have to find a prescription for 
these typical four dimensional objects. Here we 
choose the prescription of HVBM [5] or equiva- 
lently the one given by [6]. Choosing the latter 
we obtain the prescription 

1. Replace the 75-matrix by 



7^75 = g Y ^ 1 T ° r 



75 = ^potkY 1° 1 T I* ■ 



2. Compute all matrix elements in n dimen- 



3. Evaluate all Feynman integrals and phase 
space integrals in n-dimensions. 

4. Contract the Levi-Civita tensors in four di- 
mensions after the Feynman integrals and 
phase space integrals are carried out. 

This procedure requires an intensive tensorial re- 
duction of the loop-integrals as well as the phase 
space integrals. However the HVBM methods en- 
tails the presence of evanescent counter terms. 
They were calculated up to two-loop order in [7] . 
The result has the following structure 



<H 1 - •'•) T »* ( ^2" 



;/2 



-qq 



+ez, 



(1) 

1<l 



,(2),NS,- 



(11) 



where a s = a s (/i 2 )/47r is the renormalized cou- 
pling constant. Further z^q and (Q 2 / ' [i 2 ) ke are 
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process dependent. The other constants are pro- 
cess independent. The latter also holds for the 
following pieces 



7 5,NS, 
J qq 



i a s ~ 



„(2),NS, 
-qq 



^5,PS _-2 n2 / Q 2 
^qq - a s^e \-^2 



7 (2),PS 
-qq 



(12) 



(13) 



The constant Z qq NS ' + is chosen in such a way that 
the following relation holds 



A<- s (q s ,Q 2 //z 2 ) 



(14) 



where W^ q is the partonic structure function of 
the unpolarized reaction. We have to remove 
the spurious terms coming from the 75 and Levi- 
Civita prescription in the partonic structure func- 
tion by using these evanescent counter terms. 
Since they are present in the (anti-)quark sector 
we have to form the following products 

(Zl q y 2 AW qq , {Z^y 1 AW qg , 



(Z 5 qq )~ AW qq , AW, 



99 



(15) 



3. Results 

The processes which have to be calculated up 



to NNLO are 




1 — * 

q + q^i ■ 


(16) 


q + q -» 7* + .9 , 


(17) 


g + q(q) -» 7* + q(q) ■ 


(18) 


q + q->i* + g + g, 


(19) 


g + q(q) -» 7* + q(q) + g , 


(20) 


q + q^i* + q + q, 


(21) 


q(q) + q(q) -» 7* + q{q) + q(q) , 


(22) 
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Figure 5. The longitudinal asymmetry All plot- 
ted in percent in the range 2 < Q < 30 GeV at 
H = Q in LO, NLO and NNLO. 



g + g -> 7* + q + q ■ 



(23) 



including virtual corrections to lower order pro- 
cesses. Apart from evanescent counter terms 
which are characteristic of the 75 and Levi-Civita 
prescription the calculation proceeds in the same 
way as in the unpolarized case. The renormaliza- 
tion and mass factorization are carried out in the 
MS-scheme. We use the BB1 [8] parton density 
set to make our plots. The figures are constructed 
as follows 
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Because of the lack of the three-loop anomalous 
dimensions we have no Af NNL °^ Therefore in 
NNLO we use the NLO parton density set and 
the two-loop corrected running coupling constant 
a NLO F ur t;h er we set the mass factorization scale 
equal to the renormalization scale unless stated 
otherwise. Since we also compare with the unpo- 
larized Drell-Yan process we choose the following 
parton density sets given by [9]. Since they have 
an approximate NNLO set we plot the figures ac- 
cording to 
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Our findings are that the quark-anti-quark chan- 
nel dominates the proton-proton reaction in po- 
larized as well as well in unpolarized physics. The 
invariant mass distribution is shown in Fig. 3 in 
LO, NLO and NNLO. The K-factors defined by 



LO 



A(7 LO 



(24) 



are roughly the same. For Q = 7 GeV at v S = 
200 GeV we get a minimum. Here the values are 
AK NLO = 1.2 and AK NNLO = 1.3 respectively. 
In Fig. 4 we have plotted for Q = 5 GeV the 
mass factorization scale dependence 



N 



Acr(^o) ' 



(25) 



for no = Q and 0.4 fi < [i < 2 fi . We see a clear 
improvement while going from LO to NLO and 
surprisingly also from NLO to NNLO. In Fig. 5 
we have shown the longitudinal asymmetry 



A 



LL 



(26) 



in LO, NLO and NNLO. To compare the different 
polarized parton densities we have also plotted in 
Fig. 6 the longitudinal asymmetry for the BB2 [8] 
and the GRSV standard (SS) and valence (VS) 
parametrizations [10]. The figure reveals that the 
largest gluon (BB1) leads to smallest asymme- 
try. On the other hand the smallest gluon (VS) 
gives the largest asymmetry. Notice that the lat- 
ter paramctrization produces the smallest asym- 
metry in the differential p^-distribution. 
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